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Impulse  formalism  for  atom-diatom  collisions 
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(Received  15  June  1990) 

An  exact  lormulation  of  the  impulse  approach  (IA),  or  quantum-mechanical  spectator  model,  is 
applied  to  atom-diatom  collisions.  Results  are  compared  with  previous  work  on  the  IA,  which  has 
always  involved  the  peaking  approximation  (PA).  The  PA  is  seen  to  overestimate  (underestimate) 
differential  cross  sections  for  processes  involving  projectile  atom  energy  loss  (gain).  The  internal 
consistency  of  the  IA  is  explored  by  subjecting  it  to  semidetailed  balancing.  For  small  scattering  an¬ 
gles  the  IA  is  seen  to  be  an  inadequate  theory,  probably  due  to  the  neglect  of  double-  and  higher- 
collision  terms  in  the  multiple-collision  expansion  of  the  three-body  T  matrix.  For  large  scattering 
angles,  where  the  IA  does  appear  to  describe  the  scattering  process  accurately,  the  exact  calculation 
is  shown  to  give  the  same  results  as  when  only  the  energy-conserving  on-the-energy-shell  two-body 
processes  are  considered.  An  accurate  approximation  method  is  also  developed  for  rapid  computa¬ 
tion  of  inelastic  differential  cross  sections.  Finally,  the  calculated  results  are  compared  with  the  ex¬ 
perimental  measurements,  and  the  need  to  explore  two-body  potentials  more  complicated  than  the 
hard-core  potential  is  pointed  out. 


I.  INTRODUCTION 

The  basic  idea  of  the  impulse  approach  to  analyzing 
collisions  was  proposed  by  Chew1  in  his  study  of  inelastic 
scattering  of  high-energy  neutrons  by  deuterons.  When 
the  time  duration  of  the  collision  is  much  shorter  than 
the  period  of  characteristic  motion  of  the  bound  particle, 
the  components  of  the  bound  particle  may  be  considered 
to  act  independently  of  each  other,  with  position  or 
momentum  amplitudes  given  by  the  bound-state  wave 
function.  The  binding  potential  only  serves  to  create  a 
relative  momentum  distribution  of  target  particles.  The 
three-body  scattering  amplitude  thus  reduces  to  an  ap¬ 
propriately  weighted  sum  of  two-body  scattering  ampli¬ 
tudes.  The  constituent  not  taking  part  in  the  scattering 
process  is  called  the  spectator.  The  impulse  approach 
has  therefore  also  become  known  as  the  quantum- 
mechanical  spectator  model.  The  scattered  state  is  ob¬ 
tained  by  combining  the  effects  of  the  two  two-body  col¬ 
lisions  with  proper  phase  factors.  The  probability  ampli¬ 
tude  of  a  given  final  state  is  obtained  by  projecting  the 
scattered  state  onto  it.  It  should  be  emphasized  that  the 
perturbation  of  the  initial  state  is  not  assumed  small;  the 
only  assumption  is  that  the  two  scattering  centers  act 
suddenly,  and  therefore  independently.  Formal  develop¬ 
ment  of  Chew's  work  has  been  provided  by  Chew  and 
Wick,2  Ashkin  and  Wick,'  and  Chew  and  Goldberger.4 
The  impulse  approach  has  been  discussed  in  texts  by  Gol- 
berger  and  Watson'  and  Rodberg  and  Thaler,'’  the  latter 
giving  a  very  readable  account.  The  impulse  approach 
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(IA)  has  been  applied  to  a  number  of  situations  in  scatter¬ 
ing  theory.  Coleman1  has  developed  the  IA  with  applica¬ 
tions  to  electron-atom  collisions  and  Korsch  et  al. s  have 
developed  the  IA  with  applications  to  electron-molecule 
collisions.  It  has  also  been  applied  extensively  in  nuclear 
physics.5 

The  IA  was  first  applied  to  atom-diatom  collisions  by 
Bogang  and  by  Eckelt,  Korsch,  and  Philipp  (EKP).10  The 
latter  authors  in  a  series  of  papers"  "  have  applied  the 
IA  to  a  number  of  situations  involving  atom-diatom  col¬ 
lisions.  The  IA  has  also  been  applied  to  atom-diatom  col¬ 
lisions  by  Beard  and  Micha.14  In  the  I  A,  the  T  matrix  is 
expanded  in  a  multiple-collision  series111  and  only  the 
terms  representing  single  collisions  are  retained  —  two 
terms  for  an  atom-diatom  collision  and  so  on.  In  addi¬ 
tion,  while  not  a  part  of  the  basic  IA  formalism,  the 
atom-diatom  calculations  cited  above  use  another  ap¬ 
proximation  as  well,  namely  the  peaking  approximation 
(PA),  to  obtain  results  for  comparison  with  experiments. 
In  the  PA,  (i)  the  t  matrix  representing  the  two-body  pro¬ 
cess  is  evaluated  at  a  particular  value  of  the  spectator 
atom  momentum,  and  (ii)  the  slow  variation  of  the  t- 
matrix  elements  over  a  region  of  appreciable  overlap  of 
the  initial  and  the  final  states  is  used  to  justify  the  separa¬ 
tion  of  the  three-body  T-matrix  elements  into  two 
parts — one  containing  information  about  dynamics  (the 
two-body  t  matrix*  and  the  other  containing  information 
about  the  target  (the  form  factor).  7’1'1  will  be  used  to 
denote  the  three-body  transition  matrix  while  the  two- 
body  transition  matrix  will  be  denoted  by  The  PA 
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gives  a  very  nice  physical  picture  of  the  impulse  scatter¬ 
ing  process.  However,  it  was  recently  shown  by  the 
present  authors15  that,  except  for  elastic  scattering,  the 
PA  results  are  substantially  different  from  the  exact  IA 
calculation.  We  note  here  in  passing  that  the  PA  was 
also  found  to  be  inadequate  for  the  excitation  of  atomic 
hydrogen  by  1-500-keV  protons.16 

In  this  paper  we  give,  in  detail,  an  exact  formulation  of 
the  impulse  calculation  without  invoking  the  PA.  Also, 
as  in  past  studies,  the  target  molecule  has  been  taken  to 
be  in  a  low  state  of  internal  excitation,  and  a  large 
translational  energy  of  relative  motion  is  assumed  to  jus¬ 
tify  the  use  of  'he  IA  and  the  PA.  It  was  pointed  out  by 
Sharma17  that  the  IA  can  also  be  used  in  situations  in¬ 
volving  large  internal  excitation  and  small  or  moderate 
relative  translational  energies.  Since  our  formalism  does 
not  invoke  any  special  properties  of  the  T  matrix  or  the 
wave  functions,  the  resulting  equations  are  applicable  for 
arbitrary  internal  molecular  and  relative  kinetic  energies. 
Section  II  formulates  the  problem,  where  we  follow  the 
methodology  of  Eckelt  et  al. 10  Section  III  derives  the 
peaking  approximation  to  relate  the  present  work  to  the 
previous  work  and  to  set  the  stage  for  the  exact  calcula¬ 
tion,  which  is  derived  in  Sec.  IV.  Section  V  compares  the 
exact  results  with  the  PA  results.  Section  VI  checks  t  le 
internal  consistency  of  the  results  obtained  from  the  IA, 
using  semidetailed  balancing.  Section  VII  gives  an  ap¬ 
proximation  for  the  scattering  amplitude,  which  is  easily 
computed,  yet  accurate.  Section  VIII  compares  the  cal¬ 
culated  results  with  the  experimental  measurements. 
Concluding  remarks  are  given  in  Sec.  IX. 

II.  FORMULATION 

The  differential  cross  section,  in  the  center-of-mass  sys¬ 
tem,  for  a  rotational-vibrational  transition  is  given  by10 


In  a  coordinate  system  with  the  origin  coincident  with 
the  c.m.  of  the  three  particles,  only  two  momenta  are  in¬ 
dependent,  and  any  p  or  q  may  be  expressed  in  terms  of 
two  other  p’s  and/or  q’s.  The  vibrational-rotational  en¬ 
ergy  of  the  molecule  12  will  be  expressed  through  q„ 
which  will  later  be  written  in  terms  of  q,  (q2)  and  p,  (p2) 
when  atom  1  (2)  is  the  spectator  atom.  The  spectator- 
atom  momentum  remains  constant  during  the  collision. 
At  the  end  of  the  collision  qj  is  resynthesized  from  q't  and 
the  spectator  momentum  p'  =p,,  s  =  1,2.  The  probabili¬ 
ty  amplitude  of  the  final  state  is  the  projection  of  the  state 
<  q3 1  on  the  final  state. 

It  has  been  pointed  out  that  the  validity  of  the  impulse 
approach  requires  that  the  time  duration  of  the  collision 
be  much  shorter  than  the  period  of  characteristic  motion 
of  the  target.  Previous  work  on  inelastic  atom-diatom 
scattering  using  the  IA  has  assumed  a  hard-sphere  or 
hard-core  potential  between  the  incident  atom  and  either 
of  the  two  atoms  constituting  the  target.  This  potential 
satisfies  the  IA  criterion  for  any  internal  or  relative 
translational  energy  of  motion.  Another  reason  for  using 
this  model  potential  is  that  the  required  half-on-shell,  as 
well  as  off-the-shell,  two-body  /-matrix  elements  are  given 
in  literature  in  a  closed  form.  In  this  work,  we  will  also 
use  this  potential.  This  will  facilitate  the  comparison  of 
our  results  with  those  of  earlier  workers.  In  addition, 
since  our  calculation  is  exact,  we  will  be  able  to  learn 
about  the  impulse  approach  itself.  We  are,  however, 
aware  that  an  infinitely  steep  potential  overestimates 
highly  nonresonant  transitions.  In  subsequent  work,  we 
plan  to  explore  this  point  by  using  less  steep  two-body 
potentials. 

The  potential  energy  of  the  three-body  system  V  can  be 
written  as 

V  =  K,M(x,y,a!+  ,  (2) 


do 

dCl 


(vjp} 


v'j'p'i',6) 


P  3 

(2 ;  +  !)/>  3 


*22 


2 77 

ft 


2  2 
pi(d)\\T(E  +/0)|<£,>  , 
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where  |$?)  =  |uj,m,p3>  denotes  the  initial  state,  the 
final-state  variables  are  denoted  by  primes,  and  p,  and  p', 
are  the  momenta  of  the  incident  particle  before  and  after 
the  collision.  p3  is  also  the  momentum  of  the  incident 
particle  3  with  respect  to  the  center  of  mass  (c.m.)  of  the 
molecule  12.  pa  denotes  the  momentum  of  particle  a 
with  respect  to  the  c.m.  of  be,  and  qa  is  the  relative 
momentum  of  b  and  c.  This  set  of  momenta  is  called  the 
Jacobi  momenta.  denotes  the  reduced  mass  of  the  sys¬ 
tem  (a, be).  6  is  the  scattering  angle,  i.e.,  the  angle  be¬ 
tween  p,  and  p',.  The  summation  over  m  and  m'  has  re¬ 
moved  the  dependence  of  the  differential  cross  section  on 
the  azimuthal  angle.  TIE  +/0)  is  the  three-body  matrix 
element  with  total  initial  (or  final)  energy  E.  fiuh  is  the  re¬ 
duced  mass  of  the  two  particles  a  and  b ,  and 
M  -  m  j 


where  is  the  interaction  between  incident  atom  3 
and  the  diatom  M.  x  is  the  distance  of  the  incident  atom 
3  from  the  c.m.  of  the  molecule,  y  is  the  internuclear  dis¬ 
tance  of  the  diatom,  and  a  is  the  angle  between  x  and  y. 

is  the  potential  energy  of  the  diatom  12.  To  give  a 
description  of  the  three-body  scattering  process  in  terms 
of  two-body  scattering  and  a  spectator  atom  requires  that 
the  interaction  potential  be  written  as 

KVVf  =  )-»-  K,(^2  >  •  <3) 

where  Va  =  V h(.  and  y2  (J'i  )  is  the  distance  between  atom 
1  (2)  and  atom  3, 

J*i  1*:  +  [C|  (c;)y]2  +  2xyC|  (c2)cosa|1/;:  ,  (4) 

with  c,  (c2)sm2  (m  |  )/(m  | +m2 );  the  minus  sign  ap¬ 
plies  to  subscript  1  and  the  plus  sign  applies  to  subscript 
2.  The  form  for  the  atom-molecule  potential  used  in  Eq. 
(3)  does  not  give  a  good  fit  to  the  more  common  form  of 
this  potential  in  terms  of  Legendre  polynomials.  The  two 
potentials  are  applicable  to  different  models  of  atom- 
diatom  collisions.  It  is  desirable  to  establish  a  connection 
between  the  two  models,  and  this  will  be  attempted  in  a 
future  work. 

Following  EKP,10  we  write  the  three-body  /  matrix  in 
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the  multiple-collision  series 

7’  =  7-'"-f-7’'2,4-7'n,G3r'2,+  ra)Gj7’"'+  •  •  •  ,  (5) 

where  the  three-body  transition  matrix  T's\  s  =  l,2, 
denotes  the  collision  of  atoms  3  and  t,  t  —2, 1,  s  being  the 
spectator  atom  and  G3  =  (z  —  H0  —  TjF1  being  the  prop¬ 
agator  or  the  Green’s  function  (H0  being  the  total 
kinetic-energy  operator).  Equation  (3)  has  a  simple  phys¬ 
ical  interpretation.  The  first  two  terms  denote  the  col¬ 
lision  of  atoms  32  and  31,  respectively.  The  next  two 
terms  denote  double  collision  ij,  where  3  collides  first 
with  j  and  then  with  i,  and  so  on.  It  should  be  stressed 
that  this  is  not  an  expansion  in  the  Born  series,  and  each 
of  the  terms  on  the  right-hand  side  of  Eq.  (5)  represents  a 
three-body  T  matrix.  The  1A  keeps  only  the  first  two 
terms. 

To  reduce  the  three-body  matrix  T'"  into  the  two-body 
matrix  t"\  one  writes 

r'"=K1n'1',  (6) 

where  the  Mdller  operator  H‘ "  is  defined  by 

n"'=H-(z— K,)  1 K,  ,  (7) 

Vx  =  V2i,  the  potential  between  atoms  2  and  3.  The 
Mdller  operator  in  Eq.  (7)  is  still  a  three-body  operator. 
We  now  invoke  the  spirit  of  the  IA,  i.e.,  the  time  duration 
of  the  collisions  is  much  shorter  than  periods  of  charac¬ 
teristic  motion.  This  allows  one  to  assume  that  the  inter- 
nuclear  distance  of  the  diatom  stays  constant  during  the 
collision  and  to  replace  the  operator  V},  the  potential  en¬ 
ergy  between  bound  particles  1  and  2,  by  a  number. 
Writing  z  — K,=zIA,  we  get  the  two-particle  Mdller 
operator  O,  describing  the  collision  between  particles  3 
and  2, 

n,  =  l+(z,A-H0— K.r'K,  .  (8) 

In  Eq.  (7),  z  =£Vr +£rei +<0>  where  £VR  is  the 
vibrational-rotational  energy  of  the  diatomic  molecule, 
and  £rd  ~p\  /(2/u3)  is  the  relative  atom-diatom  transla¬ 
tional  energy.  Writing  £VR  —  Vj  =<i\ /(2/uI2),  where  qy 
and  /r12  are  the  relative  momentum  and  reduced  mass  of 
atoms  1  and  2,  respectively,  we  get 

zia  =pU^fiy)  +  qj  A2fin)  +  i0  .  (9) 

zIA  is  then  the  kinetic  energy  of  the  three  particles  in  the 
center-of-mass  system.  ziA  are  the  eigenvalues  of  H0, 
eigenfunctions  being  the  free-particle  states  j <b0)  de¬ 
scribed  in  the  next  paragraph.  We  have  now  translated 
into  quantitative  language  the  assertion  of  the  IA  that  the 
role  of  the  diatom  potential  is  to  generate  a  relative 
momentum  distribution  of  its  two  constituent  atoms.  We 
should  point  out  that  there  are  other  choices  for  z!A  that 
have  been  used  in  literature.14,18 

The  transition  matrix  elements  between  eigenstates  iy 
of  //„  +  Vy  can  be  written  as 


where  the  summation  is  over  the  initial  and  final  free- 
particle  states  4>0  and  4>'0,  and  <t>y,  as  pointed  out  earlier, 
denotes  |p,)  —  lp3)® \vjnt ),  etc.  The  above  expression 
can  be  simplified  by  noting  that  |<^0)  =  |q3)»  Ip, )  pro¬ 
vides  a  complete  set  for  the  free-particle  momentum 
states, 

2  I<£oX<40I<M  =  /  lq,p,Mq.,p.'l<*}>^pVq3 

=  /!q, p,><£(q,)dq,  ,  (11) 

with  the  molecular  wave  function  in  the  momentum  rep¬ 
resentation  given  by 

(q3|^,)=^(q,) 

=  (2 rr)  i/2f  drexp(/q,T)<r|iym  > 

=  (2/7T ),/2i%(q,)YjJt})  .  (12) 

Yjm  are  the  spherical  harmonics  and  Ivj  is  defined  by 

lu{<l})=  f* dr  rlX>jlrVjl<hr)  -  <13) 

X,.j(r)  being  the  radial  part  of  the  vibration-rotation  wave 
function  and  j;  representing  the  spherical  Bessel  function 
of  order  j.  Similarly, 

=  j'«i*'(q3)t/q;,<q',p',|  ,  (14) 

where  <£'(q',)  is  the  final  molecular  state.  Using  Eqs.  (11) 

and  ( 14)  we  can  write  Eq.  ( 10)  as 

(d>\\T"'\tl>y)  =  f  <b*  (q',)<q',p',U''1q,p,> 

X<Mq,)dq,dq’,  .  (15) 

To  calculate  the  matrix  elements  of  the  two-body  tran¬ 
sition  operator  / ’ ' it  is  pertinent  to  recall  that  the  spec¬ 
tator  Jacobi  momentum  p,  stays  constant  during  the  col¬ 
lision,  and  momenta  p,  and  q,  must  be  expressed  in  terms 


of  ps  and  q,,  s  =  1,2.  qs  is  changed  to  q'  during  the  col¬ 
lision  and,  after  the  collision,  inverse  transformation  to  p', 
and  qj  is  carried  out.  These  transformations  are  effected 
by  the  relations 

qj=-|Pi+["»|/<'”i+"»2>]Pjl 

,  (16) 

q,=  |p,  +  [m:/(m| +m2)]p,|  , 

(17) 

q,  =  [m,/(m2+m,)]p1+p,  , 

(18) 

and. 

q2=  -f  [mj/fm,  +m,)]p2  +  p,| 

(19) 

Equations  (16)  and  (18)  are  used  when  atom  1  is  the  spec¬ 
tator,  while  Eqs.  (17)  and  (19)  are  used  when  atom  2  is  the 
spectator.  Expressions  for  q[  are  obtained  from  Eqs. 
( 1 6)— ( 1 9)  by  replacing  p*  by  p^.  In  a  three-atom  system. 
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with  the  origin  at  the  center  of  mass,  we  have 
Pi+p2  +  p3  =  0,  and  only  two  momenta  are  independent. 
This  choice  of  momenta  greatly  simplifies  the  solution. 

Noting  that  the  state  |p3q3  >  is  also  the  state  |psq^ ), 
s  =  1,2  [Eqs.  (16)— (19)],  we  are  in  a  position  to  write 

<  qjp'jl' iJ)lq3P3  >  =  <  p^Ips  >  <  qi  > 

=  8(p;-ps)<q.;u(s,lqs>  ,  (20) 

where  the  basic  assumption  of  the  IA  (namely  that  the 
spectator  momentum  remains  constant  during  the  col¬ 
lision)  is  reflected  by  the  Dirac  8  function.  We  can  now 


rewrite  Eq.  (15)  as 

<  <f>  3 1 TL 1 ’  |  <j>  3  >  =  f  dq'}8(  p;  -  ps  )d  q}4>* ’( ) 

X<q;|f,J'|q5>^(q3)  .  (21) 

Changing  the  integration  variable  from  q3  to  p',  analo¬ 
gous  to  the  transformation  carried  out  in  Eqs.  (18)  and 
(19),  we  get 

<f,|r‘5V,>=  f  <fq30*  (q3Xq' |/lu|q(  > <6( q3 )  .  (22) 

Going  back  to  Eq.  (1),  the  expression  for  the 
differential  scattering  cross  section  can  be  written  as 


da 

da 


(vjp}~+v'j'p'},6)  = 


P 1 


P\(2j  + 1 ) 


22 


2  77 

T 


2 

1 

\ 

p-i  U) 

\ ) 

v  =  1.2 

h) 

(23) 


At  this  point  we  note  a  particular  feature  of  the  IA.  It 
has  been  pointed  out  that  qs  changes  to  q's  during  the  col¬ 
lision  between  atoms  3  and  t.  This  change  of  the  magni¬ 
tude  of  the  relative  momentum  during  scattering  from  a 
central  field  does  not  conserve  energy  during  the  two- 
body  collisions.  After  the  collision,  only  those  combina¬ 
tions  of  p\  and  internal  molecular  energy  (v'j')  that  con¬ 
serve  energy  are  permitted.  This  nonconservation  of  en¬ 
ergy  during  the  two-body  collision  should  pose  no  prob¬ 
lem,  since  the  uncertainty  in  energy,  A E,  and  time  dura¬ 
tion  of  the  collision,  At,  are  connected  by  the  uncertainty 
relation  AE^^/Af.  In  an  impulse  collision  with  a  hard¬ 
core  potential.  At— *0.  Therefore  there  is  a  large  uncer¬ 
tainty  in  energy.  This  was  the  reason,  as  we  pointed  out 
earlier,  that  the  hard-core  potential  probably  overesti¬ 
mates  highly  nonresonant  transitions.  Energy  associated 
with  qs  is  physical  in  the  sense  that  zIA 
=  ql/(2p,u  )+p?/(2fix ),  the  kinetic  energy  of  the  three- 
body  system  in  the  c.m.  frame,  whereas  the  correspond¬ 
ing  energy  associated  with  q'^q,  is  not  physical.  For 
this  reason,  the  two-body  f-matrix  elements 
f lsl(q5,q';Z|A )  are  called  half-on-the-energy-shell  matrix 
elements.  The  use  of  momentum  q',  which  does  not  con¬ 
serve  energy,  has,  as  we  have  shown  elsewhere,18  interest¬ 
ing  consequences  for  time-reversal  (TR)  symmetry.  It 
should  be  pointed  out  that  half-on-shell  and  off-shell  ma¬ 
trix  elements  have  also  been  interpreted  as  incomplete 
collisions. 19  The  two-body  t  matrix  discussed  here  does 
not  describe  the  complete  scattering  process;  it  is  there¬ 
fore  consistent  with  this  interpretation. 

III.  THE  PEAKING  APPROXIMATION 

The  peaking  approximation  (PA)  is  based  on  the  idea 
that  the  product  of  the  initial  and  final  state  in  Eq.  (22)  is 
nonzero  in  a  finite  region  of  the  spectator  momentum  and 
the  t  matrix  is  slowly  varying  in  this  region.  Thus  we  can 
evaluate  the  t1"  matrix  element  at  some  particular  value 
of  the  momenta, 


1  PA 


=  [<q;l',s,lq5>] 


p,=p, 


(24) 


where  q3  and  p3  are  related  by  Eq.  (16)  or  (17). 

The  spectator  momentum  at  which  the  t  matrix  is  eval¬ 
uated  is  arrived  at  in  the  literature  by  the  following  argu¬ 
ment:  Since  the  momentum  distributions  at  which  initial- 
and  final-state  wave  functions,  <43(q3)  and  rf>3(q3),  are 
nonzero  only  in  the  vicinity  of 


q3=0— p,= 


m. 


(m  |  +m2 ) 


P?  . 


(25) 


and 


q3=0— *p,  =  — 


m. 


(m ,  +  m-, ) 


P* 


(26) 


the  product  of  the  wave  functions  is  substantially 
different  from  zero  only  in  the  vicinity  of  their  midpoint. 


P?=- 


m. 


2(m,+m-,) 


<p,  +  p',)  . 


(27) 


This  argument  holds  strictly  for  the  v  =  0  and  the  j—  0 
state  of  the  molecule.  If  there  is  a  large  amount  of  energy 
in  the  vibration  and/or  rotation  of  the  molecule,  the 
wave  function  is  no  longer  substantially  different  from 
zero  only  in  the  neighborhood  of  q3  =  0  (or  q3=0),  and 
the  PA  may  have  to  be  evaluated  at  some  other  value  of 
the  spectator  momentum. 

Proceeding  further,  we  can  now  rewrite  Eq.  (22)  as 

(<J3|7',,,|<J3>  =  f'PA  J </q34*'(q3)<Mq3)  .  (28) 


The  integral  in  Eq.  (28)  is  an  overlap  integral,  describing 
the  amplitude  for  the  instantaneous  transition  from  the 
initial  to  the  final  state,  due  to  impulsive  momentum  in¬ 
crement. 


q.i— q.i=q.i+(-i)' 


m. 


(29) 
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q  being  the  momentum  transferred  during  the  collision,  Taking  initial-  and  final-state  wave  functions  as 
q  =  p',— pv  This  integral  is  cast  in  familiar  form  and,  at 

~  °  U/(  r)  =  v  { r)Y  (?) 

the  same  time,  its  evaluation  is  simplified  by  a  transfer-  v  AvJ 

mation  to  coordinate  space,  leading  to  ancj 

<<£3|7’lsl|<£,>  =fpA  f  dr'P*  (r)exp(  —  /a(q-r)'P(r)  ,  (30) 


where 


a  =(  —  1  )s 


(ff!1+m2) 


^‘(T)  =  x,r(r)Yrm.(?)  ,  (33) 

(3D  respectively,  and  using  the  expansion  of  the  plane  wave  in 
spherical  harmonics,  we  get 


<^|rs,|*3>=rft4ir2±/A  fAr  C(jkj';00)C(j\j';m,[i)f,l,y  vji{q)YZll{<} )  ,  (34) 

a  4iry 

where  p  =  m'  —  m,  the  Cs  are  the  Clebsch-Gordan  coefficients,20  J  =  (2y  + 1 ),  etc.,  ±  means  that  the  plus  sign  is  taken 
for  s  =  1  and  the  minus  sign  for  s  —  2,  and 

=  fjdr  r2X^r)jK(\as\rq)xvJ(>  .  (35) 

The  differentia]  cross  section  is  obtained  by  averaging  over  the  initial  azimuthal  quantum  numbers  and  summing  over 
the  final  ones,  giving 

~(vjp^v'j'p\,e)=  If  ^ 'iT  ' kCHjkj';00) 'p1;  +(  -  1  )A7^1.  ,M/& I2  -  (36) 


an  expression  available  in  the  literature.10  The  cross  sec¬ 
tion  is  factored  into  two  parts,  the  /-matrix  part  deter¬ 
mining  the  momentum  transfer  and  the  form-factor  part 
determining  the  probability  amplitude  for  the  given 
momentum  transferred.  For  homonuclear  diatomic  mol¬ 
ecules,  the  contributions  from  the  two  target  atoms  are 
equal  in  magnitude,  and  the  phase  factors  are  such  that 
only  Ay=even  transitions  are  allowed.  The  total  cross 
section  is  obtained  from  the  differential  cross  section  by 
the  relation 

,  .  , ,  ,  2tt  cp 3 da  j 

<r(vjp}->v  j  p})=- - -J  —  qdq  .  (37) 

PiPi  'Pi-Pj1  ““ 


IV.  BEYOND  THE  PEAKING  APPROaIMATION 

It  is  useful  to  recall  that  the  two-body  /  matrix  depends 
upon  the  momenta  qs  and  q’t  and  the  angle  between  these 
momenta.10,14,21  In  terms  of  other  momenta,  this  means 
that  the  matrix  element  depends  upon  incident  momen¬ 
tum  p3,  internal  molecular  momentum  q,,  and  momen¬ 
tum  transferred  during  the  collision,  q.  Of  these,  only  p3 
and  q  are  observables.  To  proceed  beyond  the  peaking 
approximation,  the  variation  of  the  two-body  t  matrix 
must  be  properly  taken  into  account  over  the  range  of  in¬ 
tegration  in  q,  in  Eq.  (22).  This  is  nontrivial,  but  can  be 
accomplished  as  follows.  We  first  note  that  the  final-state 
wave  function  in  the  momentum  representation  is  ob¬ 
tained  from  the  final-state  wave  function  in  the  coordi¬ 
nate  representation  by  the  Fourier  transform, 


<b*  ( q\ )  =  (  277- )  ,/2  f  d rexp(  —  i'q,-r ) 


where  we  have  written  q^qj  +  c^q,  and  as  has  been 
defined  earlier  [Eq.  (31)].  We  now  define  a  function 
by  the  relation 

ipui(r)  =  l2ir)~i/2  f  d q,exp(  —  / q, •  r ) ( q'  |/ul|q,  ><i(q,) 


to  obtain 

)=  f  dril>*  (r)exp(  —  /atqT)t/>u'(r)  .  (40) 

Now  we  expand  the  two-body  /  matrix  in  a  spherical  har¬ 
monics  expansion  in  the  angle  as 

<qj*  2)  *!.!♦/( 9.1  >p.i>q i i  >  (4i) 

L,M 

with  given  by 

ftw(03-P3.q>=/ dQi  )(q*  I f ' 'v  1  f q.t )  •  (42) 

Substituting  this  expansion  into  Eq.  (39),  and  introducing 
the  plane-wave  expansion, 

exp(-/q,T)=47i-2''  'k'YlJ?)YlJq,)j^q}r)  ,  (43) 


we  get 


L.M.K./i  4 77 A. 


XC(jLk;00)C(jLk;mM)  ,  (44) 


where 


Xexp(  —  »atq-r)t^*  (r)  , 


K‘l va < r  >  =  —  /  X dq>q \j K\ ( q  j r  )/[^  I ( q , , p  „ q  )/,., ( q , )  (45 ) 


(38) 
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and  /„y(^3 )  is  defined  by  Eq.  (13). 

Equation  (40)  has  the  form  of  an  overlap  integral,  if  < P" 
is  viewed  as  the  effective  initial  state  that  incorporates  the 
modification  to  the  original  wave  function  produced  by  a 
nonconstant  t  matrix.  In  the  peaking  approximation,  the 
q3  variation  of  <q^|tlsl|q, )  is  ignored,  and  Eq.  (39)  is  ap¬ 
proximated  by  evaluating  this  matrix  element  at  some 
specific,  fixed  value  of  q3,  leading  immediately  to  the 
simplified  form 

^ls)(r)=f p5A^(r)  ,  (46) 

and  Eq.  (40)  then  reduces  to  Eq.  (28)  of  the  preceding  sec¬ 
tion.  In  terms  of  the  spherical  harmonics  expansion  for 
flsl,  Eq.  (41),  it  is  as  if /lsl  is  isotropic,  i.e.,  only  the  L  =0, 
M  =  0  term  is  nonzero,  and  also  independent  of  the  mag¬ 
nitude  of q}.  Then, 


<q;U'sH>- 


"  1  00' 00  =  «  PA 


and  Eqs.  (44)  and  (45)  lead  to  Eq.  (46).  The  important 
point  about  the  modification  brought  about  by  the  varia¬ 
tion  of  the  matrix  element  of  t'"  over  qj  is  that  various 
terms  now  comprising  the  modified  wave  function  inter¬ 
fere  to  produce  the  final  scattering  amplitude.  So,  even  if 
the  variation  of  f'sl  over  q,  is  not  significant  and  the  non¬ 
isotropic  modifying  terms  are  small,  their  influence  on 
the  final  result  can  be  substantial.  Ignoring  this  is  the 
major  problem  with  the  PA. 

To  evaluate  the  integral  in  Eq.  (40),  we  substitute 
another  plane-wave  expansion, 

exp<  —  ia“lq-r) 

=  4n-^  (  ±i  fYfl}.(?)Yp}Ag  )jt>i  ia, \qr  )  ,  (48) 

ftr 

with  the  plus  sign  for  s  =  1  and  the  minus  sign  for  s  =2. 
Then, 


'l*3>=  2 

ft  K.  }-./<•  L.M 


(±i)pi'J~K> 

r 


Ypr  <  §  JlVi^xft*  q,Pi)C(jLk\OQ)C(jLk\mMp)C(kfij’\pym‘) 


where 

Ktokplq, p3)=  fdr  r2jp(  \a,  IrqWlSaWX'-j-lr)  .  (50) 

To  simplify  Eq.  (49)  for  computational  purposes,  it  is 
com  enient  to  choose  the  z  axis  of  the  coordinate  system 
along  q,  the  direction  of  momentum  transfer  vector. 
Then, 

„  1/2 

£  8)-o  •  «51) 

Using  the  properties  of  Clebsch-Gordan  coefficients,  we 
note  that,  in  Eq.  (49),  only  the/z  =  m'  term  is  nonzero  and 
M  =  (m'  —  m).  Noting  the  identity20 


C(jLk;m,M)C(kPj';m  +M,Q) 

=  2  W{jLj'P\kJ)(V)ul 

j 

XC(jJj';m,M)C(LpJ;M,0)  ,  (52) 

we  get  a  relation  between  the  total  angular-momentum 
change  during  the  collision,  J,  and  the  initial  and  final  ro¬ 
tational  quantum  numbers  j  and  j'.  We  also  note  the  to¬ 
tal  angular-momentum  change  J  is  composed  of  two 
parts:  L,  coming  from  the  two-body  t  matrix,  and  P,  the 
usual  component,  derived  from  the  momentum 
transferred  during  the  collision.  Substituting  Eq.  (52) 
into  Eq.  (49),  we  get 


<0',ir5li<£3>=  2  (±ifitj-k) 

P.KL.J 


PN&rfqyCijLk;  00) 


XC(kPj’;00)C(jJj’\mM)C(LPJ;M0)iV(jLj'P;kJ )  .  (53) 

The  expression  for  the  differential  cross  section  is  now  obtained  by  summing  the  absolute  square  of  the  collision  am¬ 
plitude  over  m  and  m'.  Since  M  =m'  —  m,  one  can  instead  sum  over  m  and  then  over  M.  Summation  over  m,  for  fixed 
M,  leads  to 


=  L. 

J 


where  the  second  factor  on  the  left-hand  side  comes  from  the  complex  conjugate  of  Eq.  (53).  Carrying  out  the  m  and 
the  J  sums,  we  finally  obtain 


2n  t  P\J  ,J  *  >' 

-f  2  2  l^.wl  . 

n  j  j  r  m 
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where 


Fjm 


2  u){j~k+0'$ 

up.  k 


a 

4ir 


1/2 

C( jLk;00)C(kf}j'',00)C(LI&J \MO)W{ jLj'P\ U )[/V[^  +  (  -  1  fN1^ ]  . 


(56) 


r 


This  is  the  central  result  of  this  paper,  invoking  no  ap¬ 
proximations  in  evaluating  the  differential  cross  section 
in  the  impulse  approach,  Eq.  (23).  The  integrations  and 
summations  become  possible  by  separating  angular  parts 
of  all  relevant  expressions  and  repeated  use  of  angular- 
momentum  identities.  It  should  be  noted  that  no  as¬ 
sumptions  were  necessary  regarding  the  structure  or  the 
properties  of  the  /  matrix.  We  have  only  invoked  the 
completeness  of  the  spherical  harmonics  expansion,  Eq. 
(41),  in  arriving  at  this  result.  Equation  (55)  then  consti- , 
tutes  the  full  calculation  of  the  differential  cross  section 
in  the  impulse  approach  for  any  two-body  potential. 

To  see  how  Eqs.  (55)  and  (56)  reduce  to  the  PA,  we 
note  that,  in  the  PA,  t"  is  assigned  some  fixed  value  in¬ 
dependent  of  q„  i.e., 

f(qj,p.„q)— fl?’(p3,q)pA  • 


Thus  L  —M  =  0  is  the  only  nonzero  term  in  Eqs.  (41)  and 
(42),  and 


/ 


(VI 

00 


—  V7  4777 


(51  . 

PA  > 


(57) 


all  other  t\ V.  =0,  and  we  can  set  L  =0  and  M  =  0  in  Eq. 
(56).  Then, 


C(  0/3./;  00)  =  Spy  , 
C(y0A.;00)  =  5M  , 


and 


■■ 

lyOOjJ 


(58) 

(59) 

(60) 

(61) 


,7=0. 

When  Eqs.  (57)— (61 )  are  substituted  in  Eq.  (55),  we  get 


all  other  NLMjJ  =0. 


j^lvjpy^v'j'p'i,e)  = 


2tt 

a 


«7  'I'  C2(jpj';00)$\f'ry  ,7  .yjf  pa  +  (  -  1  I2 

Px  ft-  r-j 


(62) 


Noting  that  /3  has  the  same  parity  as  Ay,  we  recover  the 
PA  formula,  Eq.  (36).  We  have  written  =>  instead  of  = 
in  Eq.  (61)  to  emphasize  that  Nls)  is  based  on  the  aver¬ 
aged  isotropic  component  of  rm,  whereas  the  PA  evalu¬ 
ates  at  a  given  fixed  value  of  the  spectator  momen¬ 
tum,  and  thus  ignores  the  dependence  on  the  molecular 
momentum  q3.  One  can  take  advantage  of  the 
simplification  afforded  by  the  angular-momentum  algebra 
without  completely  sacrificing  the  dependence  of  f"1 
by  taking  just  the  isotropic  component  of  the  t  matrix, 
leading  to 

N$jj* 0  ,  (63) 

all  other  N'L'y =0.  We  again  recover  the  PA  formalism, 
with  N'wjj  replacing  v7 47 Tf"'tp^  in  Eq.  (62).  We  might 
call  this  the  isotropic  approximation.  Although  the  alge¬ 
bra  involved  becomes  the  same,  the  numerical  results 
may  be  quite  different. 


V.  COMPARISON  OF  THE  EXACT  AND  PA  RESULTS 

In  this  section,  we  compare  the  results  obtained  by  our 
exact  calculation  [Eq.  (55)]  with  those  obtained  by  using 
the  PA  [Eq.  (36)].  We  use  both  formulations  to  compute 
the  differential  and  total  cross  sections  for  the  much- 
studied  process*  14,15 

Li  *  +  N2(i\y)  —  Li  *  +  N2(u',y’)  .  (64) 


For  homonuclear  diatoms,  the  two  amplitudes  corre¬ 
sponding  to  s  =  1,2  are  equal  in  magnitude.  The  relative 
phase  of  these  amplitudes  can  be  established  by  noting 

<q',|/(  1  )lq,  >(q,)  =  <q'Jf(2)|ql  >(  —  q,)  (65) 

giving 

<66) 

and 

=  (  —  1  )LN'l~m-,p  .  (67) 

Since  L  +(3  must  have  the  same  parity  as  Ay  for  the 
Clebsch-Gordan  coefficients  in  Eq.  (56)  to  be  nonzero,  the 
last  factor  in  that  equation  simplifies  to 

+  .  (68) 

Thus,  Ay  =odd  transitions  vanish,  and  the  calculations 
for  the  Ay  =  even  processes  are  simplified. 

The  two-body  interaction  is  represented  by  a  hard-core 
potential,  for  which  the  t  matrix  is  available  in  a  closed 
form.21  Since  this  is  the  only  potential  used  in  the  previ¬ 
ous  studies,  a  direct  comparison  of  our  results  with  those 
of  previous  studies  is  possible.  The  two-body  f-matrix 
element  is 

<q' |fM'(e,  )|q,  >  =/"'(q',q,,e, )  ,  (69) 

where  e,  =  g2/2jr,u  is  the  relative  translational  energy  of 
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the  two  colliding  atoms,  t  and  u,  and  is  related  to  the 
three-body  kinetic  energy  E0  by  the  relation 

E0=p]/2ni  +  q]/2fin 

=P52/2/is+952/2/xtu  •  (70) 

The  matrix  elements  in  Eq.  (69)  are  called  half-on-energy 
shells  since  the  energy  of  the  final  state  ( q'  |  is  not  equal 
to  that  of  the  initial  state  |qs ),  which  has  the  physical  en¬ 
ergy.  This  choice  of  the  t  matrix  is  called  the  “post” 
form.10  The  other  choice  of  the  t  matrix  rlsl(q',qs,e' )  is 
called  the  “prior”  form.  Since  e'  =q^2/2filu¥=ts,  the  two 
sets  of  matrix  elements  can  differ  substantially  under 
some  circumstances.  This  difference  has  been  termed  the 
“post-prior  discrepancy.” 

The  needed  vibration-rotation  wave  functions  for  N2 
were  obtained  by  numerical  integration  of  the 
Schrodinger  equation.  The  intramolecular  potential 
function  used  for  N2  was  constructed  from  spectroscopic 
constants22  using  the  Rydberg-Klein-Rees  method.23 

The  computations  are  carried  out  as  follows.  The  ma¬ 
trix  element  f 1 1  *( q', , q , ,  e , )  is  evaluated  using  Eqs.  (33)  and 
(39)  of  Eckelt  et  a/.,10  and  the  partial  waves  are 
summed  until  the  desired  convergence  is  attained.  We 
use  the  exact  expression  for  the  angle  between  q,  and 
q'l  =q,  +q.  The  range  of  interest  for  q,  is  given  by 

q,  =  j  1  —  [m,  /(m  ,  +m2)][mi/{m2  +  m}  )]|p3 

-[ffij/(m2-l-mj)]qj  ,  (71) 

obtained  by  eliminating  p2  from  Eqs.  (16)  and  (18).  It  is 
convenient  to  choose  the  coordinate  axis  so  that  q,  the 
momentum  transferred,  is  along  the  z  axis;  p3  is  taken  in 
the  z-x  plane.  p',  =  p3  +  q,  then,  is  also  in  the  z-x  plane  at 
the  center-of-mass  scattering  angle  6.  The  range  of  q,  is 
determined  by  the  wave  function  of  the  initial  state  in  the 
momentum  space,  <6( q3 ).  The  coefficients  tLV1{q}, p 3, q ) 
are  determined  by  integrating  the  product  of 
/"’(q'l.q^e,)  and  F*M(§3)  over  the  direction  of  q3,  Eq. 
(42).  Subsequent  integration  over  the  magnitude  qy  leads 
to  Klmi,  the  radial  component  of  the  wave  function  of 
the  effective  initial  state  t i.e.,  the  initial-state  wave 
function  modified  by  the  t  matrix,  Eqs.  (44)  and  (45).  In¬ 
tegration  of  over  r,  weighted  by  the  amplitude  of 
the  momentum  transferred,  q,  and  the  final-state  wave 
function,  Eq.  (50),  gives  the  dynamical  component 

of  the  transition  matrix  between  initial  and  final  states, 
Eq.  (53). 

The  typical  number  of  points  required  to  obtain  a 
significant  sixth  figure  is  as  follows:  64  points  for  the  in¬ 
tegration  over  the  direction  of  §3  (eight  points  for  the  po¬ 
lar  angle  and  eight  points  for  the  azimuthal  angle),  120 
points  for  the  magnitude  of  qy,  and  70  points  for  the  in¬ 
tegration  over  r. 

Now  we  study  the  excitation  of  N2  from  the  ground 
state  (u  =0  and  j  =0)  during  collision  with  Li  +  .  The  ex¬ 
act  and  PA  results  for  this  process,  Eq.  (64),  at  a  relative 
translational  energy  of  1  eV,  are  shown  by  plotting  the 
differential  cross  section,  at  a  c.m.  scattering  angle  of  30°, 
as  a  function  of  energy  loss  A e  =  EreI  —  £'e(  =(p2  —  pj2 )/ 


FIG.  1.  Differential  cross  section  (cm2)  for  the  exact  and  PA 
calculations  vs  energy  loss  (eV)  for  the  c.m.  scattering  angle  of 
30°  and  for  the  collision  process  Li +  -t-Nitu  =0,j  =0) 
—►Li f  +N2(u',/),  at  an  initial  relative  translational  energy  of  1 
eV.  The  peaking  calculation  (O  for  v’=0  and  A  for  u'=l) 
gives  a  larger  cross  section  than  the  exact  calculation  (solid 
symbols). 

2/Xj,  for  v’=0  and  1  in  Fig.  1  and  v’  =  2  and  3  in  Fig.  2. 
The  same  differential  cross  section,  at  a  c.m.  scattering 
angle  of  150°,  is  plotted  as  a  function  of  energy  loss  for 
u'=0  and  2  in  Fig.  3,  and  for  v'=  1  and  3  in  Fig.  4.  The 
hard-core  radius  for  the  Li  +  -N  two-body  collision  is  tak- 
en  to  be  1.62  A  for  this  energy.  For  deexcitation  pro¬ 
cesses  (Ae  <0),  we  find,  in  Fig.  5,  that  the  PA,  except  for 
small  angles  (where  the  I A  is  of  dubious  validity19),  un¬ 
derestimates  the  exact  IA  results.  Thus,  we  conclude 
that  the  PA  overestimates  (underestimates)  the  exact  IA 
cross  section  for  the  excitation  (deexcitation)  process.  It 
is  interesting  to  note  that  both  the  PA  and  exact  calcula¬ 
tions  lead  to  the  same  general  pattern  for  the  differential 
cross  section.  The  oscillatory  structure  of  the  state-to- 
state  differential  cross  sections,  especially  the  sharp  final 
peak  calculated  for  each  v'  termed  “rotational  rainbow 
structure,”8  is  exhibited  both  by  the  PA  and  the  exact 
calculation.  This  similarity  of  patterns  is  a  reflection  of 


FIG.  2.  Differential  cross  section  (cm2)  for  the  exact  and  PA 
calculations  vs  energy  loss  (eV)  for  the  c.m.  scattering  angle  of 
30°  and  for  the  collision  process  Li  *  +N,(e  =0,y  =0) 
-»Li '  +N,(t>',/),  at  an  initial  relative  translational  energy  of  1 
eV.  The  peaking  calculation  (O  for  n'  =  2  and  A  for  r'  =  3) 
gives  a  larger  cross  section  than  the  exact  calculation  (solid 
symbols). 


43 


•  IMPULSE  FORMALISM  FOR  ATOM-DIATOM  COLLISIONS 


197 


FIG.  3.  Differential  cross  section  (cm’l  for  the  exact  and  PA 
calculations  vs  energy  loss  (eV)  for  the  c.m.  scattering  angle  of 
150°  and  for  the  collision  process  Li "  +  N,(e  =0,j  =0) 
— -Li'  +N2(i'',/),  at  an  initial  relative  translational  energy  of  1 
eV.  The  peaking  calculation  (C  for  u'=0  and  A  for  e'  =  2) 
gives  a  larger  cross  section  than  the  exact  calculation  (solid 
symbols).  Notice  the  larger  rotational  excitation  than  that  for 
30°  scattering. 

the  fact  that  small  L^0,M  =0  terms  in  Eq.  (56)  interfere 
with  the  dominant  isotropic  term.  It  has  been  pointed 
out  earlier15  that  keeping  only  the  isotropic  (L  =0)  term 
has  the  same  mathematical  structure  as  the  PA  but  with 
a  different  t  matrix.  This  interference  between  different  L 
terms  modifies  the  numerical  value  of  the  differential 
cross  section,  keeping  the  general  pattern  intact.  This 
point  is  further  discussed  in  Sec.  VII. 

Finally,  we  display  the  total  cross  section 
CTr(00=>u7')  for  the  Li++N2  collisions  at  relative 
translational  energies  of  1  eV  (Fig.  6)  and  10  eV  (Fig.  7). 
It  should  be  noticed  that  for  vibrationally  inelastic  col¬ 
lisions  the  cross  section  increases  rapidly  at  the  thresh¬ 
old,  stays  constant  as  higher  rotational  levels  are  popu¬ 
lated  with  greater  energy  loss,  and  then  falls  rapidly.  For 
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FIG.  4.  Differential  cross  section  (cnr)  for  the  exact  and  PA 
calculations  vs  energy  loss  (eV)  for  the  c.m.  scattering  angle  of 
150°  and  for  the  collision  process  Li '  +  N,(i>  =0,y  =0) 
-•Li '  +N,(r',/),  at  an  initial  relative  translational  energy  of  1 
eV.  The  peaking  calculation  (  o  for  v'  =  1  and  A  for  i>'  =  3)  gives 
a  larger  cross  section  than  the  exact  calculation  (solid  symbols). 
Again,  notice  the  larger  rotational  excitation  than  that  for  30° 
scattering. 
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FIG.  5.  Differential  cross  section  (cm:l  for  the  exact  and  PA 
calculations  as  a  function  of  the  center-of-mass  scattering  angle 
(deg)  for  the  deexcitation  process  Li '  +  N,(6,6)— Li ' 
+  N,(  3, 12)  at  a  relative  translational  energy  of  1  eV. 

1-eV  relative  translational  energy  and  final  vibrational 
levels  N2(u’  =  2  and  3),  this  fall  occurs  when  the  available 
energy  is  almost  totally  depleted.  This  result  is  consistent 
with  the  measurements  of  Loesch  and  Herschbach24  who, 
in  their  work  on  crossed  be?ms  of  Ar  and  Csl  at  collision 
energies  of  0.35-1.1  eV,  found  peaks  that  “correspond  to 
an  extremely  inelastic  ‘ballistic’  process  in  which  most  of 
the  initial  relative  translational  energy  goes  into  vibra¬ 
tional  or  rotational  excitation."  King,  Loesch,  and  Her¬ 
schbach  reported25  similar  observations  from  a  crossed- 
beam  study  of  Ar  and  CsF  at  a  relative  kinetic  energy  of 
0.66  eV.  For  N2(u'  =  0),  the  excitation  cross  section 
drops  off  dramatically  when  about  half  of  the  available 
translational  energy  goes  into  rotational  energy.  This  is 
because  when  momentum  q  is  transferred  to  the  mole¬ 
cule,  its  rotational  quantum  number  may  be  altered  by  up 
to  |a,  I bq,  where  b  is  the  internuclear  distance  of  the  dia¬ 
tom,  to  conserve  total  angular  momentum.  Further,  if 
the  molecule  is  initially  moving  slowly  or  is  at  rest  in  the 
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FIG.  6.  Total  cross  section  according  to  (he  exact  calculation 
as  a  function  of  energy  loss  for  the  collision  of  N;  with  Li  ‘ ,  at  a 
relative  translational  energy  of  1  eV.  The  molecule  is  initially  in 
the  v  =0  and  j  =0  states.  Final  vibrational  states  are  denoted  as 
follows:  ...  r'  =  0;  \,  v'=  1;  □,  t>'  =  2;  and  0.  r'  =  3.  Final  rota¬ 
tional  states,  /,  increase  sequentially  from  left  to  right,  starting 
with  /= 0  in  steps  of  2.  Notice  that  almost  all  the  translatio"',l 
energy  may  be  converted  into  internal  degrees  of  freedom. 
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FIG.  7.  Same  as  Fig.  5,  but  with  £r<;1  =  10  eV.  Several  final 
vibrational  states  are  shown,  with  symbols  as  defined  in  the 
figure.  A  precipitous  drop  in  the  excitation  cross  sections  is  to 
be  noticed. 

laboratory  frame,  it  acquires  [1/2 (m{+m2)]q2  kinetic 
energy  to  conserve  the  total  linear  momentum.  The 
physical  picture  that  emerges  is  that  if  the  momentum 
transferred  is  comparable  to  the  momentum  of  the  center 
of  mass  of  the  molecule,  we  observe  the  “ballistic” 
phenomenon.  On  the  other  hand,  when  the  momentum 
transferred  is  much  larger  than  the  c.m.  momentum  and 
the  molecule  is  excited  to  very  high  rotational  or  vibra¬ 
tional  levels,  not  all  the  initial  relative  translational  ener¬ 
gy  can  be  converted  into  internal  energy.  This  feature 
was  evident  in  the  experiment  of  Hershberger  el  al.,2*' 
who  observed  that  during  the  collisions  of  C02  gas  at 
room  temperature  with  energetic  H,  D,  and  Cl  atoms,  the 
molecules  excited  to  higher  rotational  levels  are  also  the 
ones  moving  at  greater  speeds. 


VI.  INTERNAL  CONSISTENCY 
OF  THE  CALCULATED  CROSS  SECTIONS 

The  previous  sections  have  shown  how  to  do  an  exact 
calculation  of  the  inelastic  atom-diatom  cross  sections  us¬ 
ing  the  IA  formalism.  In  this  section,  we  check  the  re¬ 
sults  for  internal  consistency  by  testing  them  for  semide- 
tailed  balancing  (SDB). 

The  Hamiltonian  is  independent  of  time  and  quadratic 
in  momenta,  and  is,  therefore,  invariant  under  time  rever¬ 
sal  (TR).  This  symmetry,  together  with  space  invariance, 
requires27  that  the  exact  scattering  amplitudes  obey  the 
SDB  relation  between  the  differential  cross  sections  for 
the  forward  and  reverse  processes,  i.e., 

p[da(i  -+f)/dSl]/p(  f)=p'[da(f-—i)/dn]/p(i)  , 

(72) 

where  p(i)  and  p( /)  are  the  densities  of  state  in  the  initial 
and  final  states;  p  and  p'  represent  the  incident  flux  densi¬ 
ties  for  the  direct  and  reverse  processes.  It  was  shown 
earlier18  that  the  differential  cross  sections  for  the 
rotational-vibrational  excitation  of  N2  upon  collision  with 
Li  *  do  not  satisfy  SDB  for  small  scattering  angles.  The 
scattering  amplitudes  for  the  forward  process  using  the 


post  model  of  the  two-body  energy  parameter  and  the  re¬ 
verse  process  using  the  prior  model  of  the  two-body  ener¬ 
gy  parameter  are  the  same  (Appendix  A),  apart  from  an 
overall  sign  factor  and  change  in  the  sign  of  the  magnetic 
quantum  numbers  ( m  — *  —  m  and  m '  — ►  —  m ' ).  We  have, 
therefore,  the  general  result  that  the  differential  cross  sec¬ 
tions  for  the  forward  and  the  reverse  process  are  related 
by  the  density-of-states  factors.  Thus,  the  failure  of  the 
post  model  to  satisfy  SDB  also  implies  the  corresponding 
failure  of  the  prior  model.  That  these  two  models  are  not 
adequate  for  describing  the  collision  process  for  all  angles 
is  shown  by  the  fact  that  the  total  scattering  cross  section 
using  these  models  for  two-body  energy  is  much  less  than 
the  geometric  shadow.  The  ‘full”  model14  of  the  two- 
body  energy  parameter  z1A  obeys  SDB  for  all  angles,  but, 
as  has  been  noted  earlier,18  suffers  from  other  difficulties, 
one  being  that  the  total  calculated  cross  section  is  far  in 
excess  of  the  geometric  shadow. 

For  large  scattering  angles,  the  differential  cross  sec¬ 
tion  was  shown18  to  be  rather  insensitive  to  the  model 
used  for  the  two-body  energy  parameter  z,A.  This  can  be 
understood  as  follows.  The  two-body  scattering  ampli¬ 
tude  from  a  hard-core  potential  for  large  scattering  an¬ 
gles,  provided  qsrc»  1,  is  constant,  independent  of  col¬ 
lision  energy,  as  long  as  the  ratio  of  the  incident  to  the 
final  momentum  stays  the  same  (Fig.  8)  (the  two-body 
scattering  angle  is  denoted  by  <f>  while  the  three-body 
scattering  angle  is  represented  by  6).  On  the  energy  shell, 
i.e.,  q2  =  q^2,  this  value28  is  rt./ 2.  The  equivalent  classical 
process  would  be  the  energy-conserving  on-the-energy- 
shefl  process,  one  with  ratio  q ,  /q'  =  1.  Using  the  relation 
q'  =  —  (  —  1  )'q  +  q,,  and  recalling  that  the  z  axis  in  this 
calculation  is  taken  along  the  momentum  transfer  vector 
q,  the  on-the-energy-shell  requirement  becomes 
—  (  —  1  V2q,,q  +  <?2  =  0,  qsz  being  the  z  component  of  q,. 
The  root  q  =0  corresponds  to  forward  scattering  with  no 
change  of  momentum,  and  in  the  classical  picture  it  cor¬ 
responds  to  no  collision.  The  root  qs:=(  —  1  leading 
to  q':  =  <?,. ,  corresponds  to  the  situation  in  which  the  z 
component  of  the  incident  momentum  changes  sign,  and 
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FIG.  8.  Absolute  value  of  the  two-body  scattering  amplitude 
divided  by  the  range  of  the  hard-core  potential,  as  a  function  of 
the  two-body  scattering  angle  <6,  for  an  incident  momentum  of 
10“  cm  '.  An  asymptote  of  0.5  is  reached  for  the  on-the- 
energy-shell  process. 
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the  x  and  y  components  remain  unaffected.  This  is  very 
much  the  description  of  a  collision  between  two  hard 
spheres  using  classical  mechanics.29  We  still  compute  the 
two-body  t  matrix  quantum  mechanically,  so  we  call  the 
on-the-energy-shell  (<j2  =  <js'2)  two-body  t  matrix  the  semi- 
classical  t  matrix.  Using  the  relation  pi;  =  —  Ae/q) 
+  (q/2),  where  A e  is  the  energy  gained  by  the  internal 
degrees  of  freedom  of  the  molecule,  and  Eqs.  (16)-(19), 
we  can  write 


qf:~  ~  [axq /2)  +  (pv_\i /axq  )  .  (73) 


Recalling  that  q}  is  the  relative  momentum  of  the  two 
atoms  in  the  diatom,  we  note  that  Eq.  (73)  restricts  the 
value  of  the  component  of  q,  along  the  momentum 
transfer  vector  q,  in  order  for  the  molecule  to  gain  energy 
Ae  by  the  energy-momentum  conservation  relations  for 
scattering  between  two  hard  spheres. 

Relation  (73)  also  results  from  quantum  mechanics  if 
q*.  is  defined  as  the  weighted  average 


9% 


f  6’ q  \:<t>dqy 
f  <b*6dq. 


(74) 


Relation  (74)  is  proved  in  Appendix  B.  Because  the 
derivation  of  Eq.  (74)  does  not  use  the  properties  of 
scattering  from  a  hard  sphere,  Eq.  (73)  is  a  more  general 
relation.  We  can  then  conclude  that  large-angle  scatter¬ 
ing  can  be  described  by  on-the-energy-shell  matrix  ele¬ 
ments  using  the  two-body  t  matrix  for  any  potential.  To 
demonstrate  the  accuracy  of  these  conclusions,  we  plot  in 
Fig.  9  the  differential  cross  section  versus  the  angle  for 
the  exact  calculation  and  the  semiclassical  calculation  us¬ 
ing  the  on-the-energy-shell  two-body  t  matrix  evaluated 
at  fixed  <?,.  =  q ?f,  given  by  Eq.  (73).  The  transverse  mo¬ 
menta  q}x  and  qM.  still  range  over  the  entire  x-y  plane. 
The  two  calculations  agree  well  at  large  angles,  showing 
that  for  these  angles  the  two-body  scattering  is 
effectively,  indeed,  an  on-the-energy-shell  process.  A 
semiclassical  calculation  restricting  q}x=qx  ~0  in  the 
two-body  t  matrix  is  also  shown  in  Fig.  9.  The  results  are 


exactly  the  same  as  from  the  calculation  in  which  qix  and 
qiv  are  not  restricted,  if  the  c.m.  scattering  angle  is  >  15°. 
This  simplification  leads  to  an  expression  for  the 
differential  cross  section  identical  to  that  given  by  the 
PA,  except  that  the  two-body  i  matrix  is  evaluated  for 
q given  by  Eq.  (73),  instead  of  q}.  =  —asq/2.  It  is  seen 
that  in  the  region  of  validity  of  the  IA,  the  two-body  t 
matrix  may  be  computed  by  the  on-the-energy-shell  semi¬ 
classical  procedure.  Further,  the  results  obtained,  at 
least  for  the  hard-core  repulsive  potential,  are  indepen¬ 
dent  of  the  model  chosen  (post  or  prior)  for  the  two-body 
energy  parameter,  z1A . 


VII.  ANALYTICAL  APPROXIMATIONS 

We  now  develop  some  analytical  approximations 
which  lead  to  a  relatively  simple  expression  for  the 
differential  cross  section,  even  for  small  scattering  angles. 
This  gives  us  further  insight  into  the  IA,  and  allows  us  to 
understand  the  similarities  of  the  exact  and  the  PA  cross 
sections  displayed  in  an  earlier  section.  In  addition,  it 
also  provides  an  effective  procedure,  valid  over  a  wide 
range  of  parameters,  for  obtaining  the  IA  cross  sections 
with  much  less  computational  effort.  The  analytical  ap¬ 
proximation  developed  in  this  section  has  an  advantage 
over  the  semiclassical  approach  of  the  last  section,  be¬ 
cause  it  is  valid  over  a  larger  range  of  scattering  angles. 
This  approximation  will  make  it  possible  to  calculate  the 
higher  terms  of  the  multiple-collision  expansion  of  the 
three-body  T  matrix,  Eq.  (5),  of  which  the  IA  constitutes 
the  single-collision  terms. 

Our  starting  point  is  the  observation  (Fig.  10)  that  the 
two-body  i  matrix  for  the  post  and  prior  models  is  a 
linear  function  of  q};.  The  corresponding  results  for  the 
full  model,  however,  cannot  be  expressed  in  so  simple  a 
form.  Recalling  that  the  value  of  q,  perpendicular  to  q 
has  no  effect  on  the  value  of  the  t  matrix,  the  z  axis  is  tak¬ 
en  along  q,  and  q\x  and  q}y  are  set  equal  to  zero.  The  ex¬ 
pression  for  the  two-body  t  matrix,  in  the  post  and  prior 
forms,  becomes 


t"  (q,)  =  (q'  t"'  q,  >  =[f  (0)  +  q}Jdt  /dq\: >()]  •  (75) 

In  mathematical  language,  the  neglect  of  q}x  and  q\y  is 
justified  by  noting  that  ( bt  /dq l„  and  (dl/dq\v\)  are 
negligible  compared  to  their  z  counterpart.  Physically, 
this  means  that  the  results  for  a  hard-core  potential  are 
insensitive  to  the  angle  between  the  plane  containing  p„ 
p’„  and  q.  and  that  containing  q,,  q', .  and  q,  so  the  two 
planes  can  be  taken  to  be  coincident.  We  can  now 
rewrite  Eq.  (22)  as 

J  A*  / 1  q -  )<Mq,==  /  (0)  J  A*  Ad q, 


FIG.  9.  Exact  and  two  semiclassical  differential  cross  sec¬ 
tions  as  a  function  of  the  ccnter-of-mass  scattering  angle,  for 
Li  FN.(O.O)  -Li '  -F  N ,<  1.6).  The  two  semiclassical  calcula¬ 
tions  become  identical  for  angles  larger  than  about  16*.  They 
also  equal  the  exact  calculation  after  about  35°. 


-Hi)/  /dqx.  ),,  J*(5*  qx.Adq}  .  )76> 

The  first  integral  in  Eq.  (76)  is  just  the  form  factor  of  the 
molecule  for  the  momentum  transfer  q. 
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F(q)=  f  <t>*  id q, 

=  f  dri*  (r)exp(  — /a,q-r)t/'(r' 

=  </k  '""'l i>.  (77) 


where  /  and  i  denote  the  molecular  final  and  initial  states. 
This  integral  was  evaluated  earlier  in  Eqs.  (28)  — (35). 
The  integral  in  the  second  term  of  Eq.  (76)  is  shown  in 
Appendix  B  to  be  qf:  of  Eq.  (73)  times  form  factor  F(q). 
The  two-body  t  matrix  now  becomes 


f  <t>*  r(q,  )idq}: 


t(0)  +  qf: 


_a t 


F(q) 


(78) 


This  is  analogous  to  the  peaking  result,  Eqs.  (28)  or  (30), 
with  t,,A  replaced  by  the  expression  in  the  large  square 
brackets.  The  remaining  algebra  for  obtaining  the 
differential  cross  section  proceeds  as  in  Eqs.  (34)  — (36). 
This  explains  why  the  shapes  of  the  PA  and  exact 
differential  cross  sections  are  so  similar. 

A  detailed  comparison  of  the  exact  and  analytical  [Eq. 
(78)]  calculations  (using  the  post  model)  is  presented  in 
Fig.  11  for  the  collision  of  N,  with  Li  * .  The  analytical 
method  of  approximating  the  differential  cross  section 
has  been  found  to  agree  with  the  exact  calculation  over  a 
wide  range  of  parameters  when  the  two-body  /  matrix  can 


FIG.  10.  (a)  Real  and  (b)  imaginary  parts  of  the  two-body 
scattering  amplitude  as  a  function  of  q,,,  the  component  of  the 
internal  momentum  along  the  momentum  transfer  vector  q  for 
post,  prior,  and  full  models  of  the  two-body  energy  parameter 

ZIA- 


FIG.  1 1.  Exact  and  analytical  [based  on  Eq.  ( 78 )]  differential 
cross  sections,  as  a  function  of  the  c.m.  scattering  angle,  for  the 
collision  of  N,  with  Li '  at  a  relative  translational  energy  of  1 
eV,  for  the  transition  (0,0  -  1,6). 


be  represented  by  the  two-term  expansion  of  Eq.  (75). 
For  the  deexcitation  processes,  the  analytical  approxima¬ 
tion  also  shows  good  agreement  with  the  exact  calcula¬ 
tion,  if  the  prior,  rather  than  the  post,  model  of  the  two- 
body  t  matrix  is  used.  The  post-prior  discrepancy  is  like¬ 
ly  to  disappear  for  smaller  angles  as  well  upon  the  in¬ 
clusion  of  the  multiple-collision  terms. 

VIII.  COMPARISON  OF  THE  COMPUTED 
AND  EXPERIMENTAL  RESULTS 

We  compare  the  results  of  our  computations  with  the 
measurements  of  Bottner,  Ross,  and  Toennies,'0  who 
crossed  a  Li f  ion  beam  with  a  N,  or  CO  beam.  Figures 
12  and  13  compare  the  total  differential  cross  section  as  a 
function  of  the  CM  scattering  angle  at  a  relative  transla- 


FIG.  12.  Measured  (solid  circles)  and  calculated  (open  cir¬ 
cles)  total  differential  cross  sections  for  the  collision  of  N,(0,2> 
with  Li  ‘ ,  at  a  relative  translational  energy  of  3.64  eV,  as  a  func¬ 
tion  of  the  center-of-mass  scattering  angle.  The  experimental 
results,  since  they  are  not  given  in  absolute  units,  are  normal¬ 
ized  to  agree  with  the  calculated  ones  at  a  scattering  angle  of 
100°.  We  note  that  a  substantial  amount  of  energy  goes  into 
higher  \ibrational  levels.  Cross  sections  for  the  final  vibrational 
levels,  summed  over  the  final  rotational  levels,  are  r'  =  0(/\), 
r'  =  It  □),  t’'  =  2(0),  t’'  =  3(V),  r'-4(*),  r’  =  5(X),  and 

t’'  =  6(  +  ). 
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FIG.  13.  Same  as  Fig.  12,  except  that  the  molecule  is  CO  and 
the  experimental  points  are  now  given  in  absolute  units,  the  nor¬ 
malization  being  provided  by  the  N,  result  of  Fig.  12.  The 
agreement  of  the  calculation  with  the  measurements  is  excel¬ 
lent. 


tional  energy  of  3.64  eV  for  N,  and  CO,  respectively. 
Measured  cross  sections  are  scaled  to  coincide  with  the 
calculated  ones  for  N,  at  a  scattering  angle  of  100°.  As 
can  be  readily  seen  from  the  figures,  the  calculated  and 
measured  cross  sections  are  in  very  good  agreement  for 
large  angles.  We  give  the  contribution  of  various  final  vi¬ 
brational  levels  to  the  total  differential  cross  section.  Fi¬ 
nally,  we  plot  (Fig.  14)  differential  cross  sections  for  49.2° 
at  a  relative  translational  energy  of  4.23  eV  as  a  function 
of  the  final  rotational  level  (y')  for  the  v'—O  and  y'=l 
levels  of  N,.  The  calculated  results  have  a  very  narrow 
distribution  and  are  peaked  at  much  higher  rotational 
levels  than  the  experimental  data.  Similar  results  were 
obtained  for  Li  *  and  CO  collisions.  These  results  are 
largely  a  reflection  of  the  hard-core  two-body  potential, 
which  gives  large  and  narrow  momentum  transfer 
(  —  1  )'(q,/2)  =  9,.  =  —  ql:.  A  more  realistic  potential  may 
give  more  reasonable  results. 


(I  in  '.’o  :io  io  .',o  no  vo 


FIG.  14.  Measured  and  calculated  differential  cross  sections 
for  the  collision  of  N:(0,2 )  with  Li  ’ ,  with  a  relative  translation¬ 
al  energy  of  7.07  eV,  at  a  center-of-mass  scattering  angle  of  49.2° 
for  final  vibrational  levels  0  and  1 .  The  peak  of  the  vibrationally 
elastic  (v'=0)  measured  cross  section  is  normalized  to  agree 
with  the  corresponding  calculated  cross  section. 


IX.  CONCLUSIONS 

A  rather  complete  description  of  the  exact  formulation 
of  the  IA  and  its  validity  is  given.  On  the  formal  side,  the 
inclusion  of  the  double  and  higher  collision  terms  in  the 
calculation,  to  be  considered  next,  may  render  the  formu¬ 
lation  self-consistent  for  smaller  scattering  angles.  It  is 
pertinent  to  point  out  here  that  the  IA  is  not  an  approxi¬ 
mation  to  fully  quantum  coupled  states  or  classical  trajec¬ 
tory  calculations.  When  the  time  duration  of  the  col¬ 
lision  is  much  smaller  than  the  period  of  the  characteris¬ 
tic  motion  of  the  diatom,  the  IA  appears  to  be  a  more  ap¬ 
propriate  theory.  In  fact,  classical  trajectory  calculations 
"so  far  have  not  been  able  to  explain  the  observed  distri¬ 
butions  in  rotational  transition  probabilities,  consistently 
yielding  distributions  that  are  too  narrow  and  peaked  at 
smaller  values  ly  «  10)  than  observed  in  the  experiments 
( j  =  30)".  11  This  conclusion  was  also  arrived  at  by  Hersh¬ 
berger  el  al.->'  More  realistic  two-body  potentials  are 
currently  being  investigated  in  an  attempt  to  use  the  IA 
to  understand  experiments  involving  large  vibrational- 
rotational  inelasticities. 

The  approach  to  vibration-rotation  energy-transfer 
processes  discussed  here  was  first  outlined  by  Landau  and 
Teller,12  who  performed  a  one-dimensional  calculation  to 
obtain  the  vibration  to  translation  rate  coefficients.  Sub¬ 
sequent  authors,  most  notably  Schwartz,  Slawsky,  and 
Herzfeld,11  have  expanded  on  the  Landau-Teller  work, 
however  still  keeping  it  a  one-dimensional  calculation. 
We  have  examined  the  basic  ideas  of  these  authors  in  the 
IA  framework,  elucidating  the  physics  involved  and  for¬ 
mulating  it  in  three  dimensions  to  demonstrate  the  im¬ 
portant  role  that  the  rotational  degrees  of  freedom  play 
in  these  processes.  Our  formulation  treats  the  internal 
momentum  of  the  diatom  due  to  vibrational  and  rotation¬ 
al  motion,  q„  on  an  equal  footing  with  the  relative 
translational  momentum,  p,.  This  makes  it  possible  for 
us  to  apply  our  calculation  to  collisions  involving  the  in¬ 
tramolecular  transfer  of  vibrational  energy  into  rotational 
energy,14  and  the  transfer  of  rotational  energy  into  vibra¬ 
tional  energy,  as  would  be  the  case  in  a  shock-tube  experi¬ 
ment.15  By  doing  an  exact  calculation  on  a  simple  poten¬ 
tial,  we  have  also  been  able  to  define  the  range  of  applica¬ 
bility  of  our  calculation  and  the  steps  needed  to  extend 
this  range. 


APPENDIX  A 

The  post  model  of  the  two-body  I  matrix  upon  time  re¬ 
versal  ITr)  changes  to  the  prior  model  and  vice  versa. 
This  is  seen  by  noting  that  incoming  (outgoing)  momenta 
become  outgoing  (incoming)  momenta  with  a  minus  sign. 
Then, 


TRt 


q,-»q,;r  = 


=  t 


=  t 


prior  ^  q* 


q,:c  = 


—  —  q,  >  • 


<?; 

^  fu 

(A  1 ) 


202 


SHARMA,  BAKSHI,  AND  SINDONI 


43 


Similarly, 

^i^prior^Qs  *<T  )  —  ,posl(  *  Q  .s  (  •  (A2) 

The  wave  functions  are  transformed  under  time  rever¬ 
sal  in  the  usual  way.27  Note,  however,  that  the  momen¬ 
tum  transfer  vector  is  invariant  under  TR, 
q  =  Pj  —  pj=  —  p,  —  (  —  pj).  Invoking  the  property  of  the 
wave  functions  in  momentum  space, 

Kh  =  ,  (A3) 

it  is  straightforward  to  see,  by  a  simple  change  of  vari¬ 
ables  in  Eq.  (22),  that  the  scattering  amplitude  for  the  for¬ 
ward  process  using  the  prior  model  of  the  two-body  ener¬ 
gy  parameter  is  related  to  that  of  the  reverse  process  us¬ 
ing  the  post  model, 

^  priori"*. P-*l’T"»',P’> 

=  .(uV'-m’, -p'  —  vj-m,-p)  . 


Using  the  space  inversion  invariance  (which  changes 
the  sign  of  linear  momenta),  we  have 

^pmi^V'-w'.-p'-*  vj-m, -p) 

=  (  —  1  V'-jA-mst{v'j'  —  m,,p'-*vj  —  m,p)  . 


Combining  Eqs.  (A4)  and  (A5),  we  get 

* poJv'j' ~ m’,p' -^vj  - m,p) 


_,),'w  J  ”A<vjm,p^vym',p'>  . 


Taking  the  absolute  square  of  both  sides  of  Eq.  (A6)  and 
summing  over  m  and  m '  leads  to  the  desired  relation  be¬ 
tween  the  forward  and  reverse  processes. 


APPENDIX  B 

We  prove  the  relationship 

c-  f  q> 

— ■ —  (Bl) 

fd>}<b,d  q, 

in  this  appendix.  We  start  by  rewriting  the  numerator  in 
Eq.  (Bl)  as 

J"  i*  <?  ,.<Mq,=  f  dttl’*  (r)exp(  — /a,q-r)§,_.tWr) 

=  ( f\e  m''qqi.\i )  .  (B2) 

This  is  proved  by  using  the  defining  relation  for  <A,(q}). 
the  molecular  wave  function  in  the  momentum  space  in 
terms  of  t/',(r),  and  the  molecular  wave  function  in  the 
coordinate  space. 


q iz <t>i ( q j >  =  ( 2 rr )  3/2 f  q i: exp[ /  ( q , +  q iyy  +  q ,.z  IJ^trWr 

=  (2rr)  3/2  j dx  dy  exp[/ (qixx  d-q)vy )]  f  dz  —i—exp(iqi2z)  i/',(r) 


Integrating  by  parts,  we  get 

f  dz  —  i~-  expdq^.z )  t/>,(r) 

J  dz 

=  —  fexp(iq}.z)  — /J^iA,(r)  dz  .  (B4) 

Writing  the  defining  relation  for  <4y(q', ),  recalling  that 
qj  =  q,4-«  q,  and  then  integrating  over  q,  we  get  Eq. 
(B2),  with 


H  =  z^-(Qlx+iily+ill)+V(rn)  , 
W 12 


where  /i)2  is  the  reduced  mass  and  rl2  is  the  internuclear 
distance  of  the  diatom,  we  can  rewrite  Eq.  (B5)  as 


[exp(  —iasqz  ),//]=  -exp(  — iasqz  )— ^ — la,q+2q^)  , 

2^I2 


<7.1.-  =' 


giving 


Then,  the  exponential  acts  as  a  translation  operator, 

exp(  —  iasqz  )/($,.  )exp(/«s^z  )=/(§,.  —  a, q)  ,  (B5) 

because  the  momentum  operator  is  conjugate  to  the 
coordinate  operator  z.  Recalling  that  the  Hamiltonian 
operator  for  the  diatom  is  (setting  ^  =  1 ) 


exp(  ~iasqz  )§,.  = - -[expi  -iasqz ),//] 

a,q 

—  exp(  —  ia,qz )-y-  .  (B8) 

The  matrix  elements  of  the  right-hand  side  of  Eq.  (B8)  be¬ 
tween  the  initial  and  final  states  give  Eq.  (Bl). 
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